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ABSTRACT 



The peculiar X-ray binary Cygnus X-3 has been observed 
on several occasions with the X/y-ray instruments on 
board INTEGRAL. We have collected data from available 
public and Galactic Plane Scan observations between De- 
cember 2002 and December 2003 and summed them to- 
gether into two broad-band spectra, representing different 
physical spectral states of the source. We have fitted the 
two spectra with models including Comptonization and 
Compton reflection, similar to those found for black-hole 
binaries at high accretion rates. 

Key words: gamma rays: observations - radiation mech- 
anisms: non-thermal - stars: individual: Cyg X-3 - X- 
rays: binaries - X-rays: general - X-rays: stars. 



1. INTRODUCTION 



Cyg X-3 is an enigmatic X-ray binary that ever since its 
discovery 38 years ago, in 1966 (Giacconi et al. 1967), 
has evaded simple classification. Its 4.8 hr orbital period 
is typical of a low-mass system, but infrared observations 
have shown the donor to be a high-mass Wolf-Rayet star 
(e.g., van Kerkwijk et al. 1992; 1996). The nature of the 
compact object is uncertain; it may be either a neutron 
star or a black hole (e.g., Schmutz et al. 1996; Ergma & 
Yungelson 1998; Stark & Saia 2003). The system, lo- 
cated at a distance of ~9 kpc (e.g., Predehl et al. 2000), is 
embedded in such a dense wind from the donor star that 
most of its intrinsic emission is strongly obscured. Still, 
its hard X-rays and y-rays can be observed unimpeded, 
making it a very good target for INTEGRAL. 

As an X-ray source, Cyg X-3 is one of the brightest in 
the Galaxy. It displays states with high and low level of 
emission at <10 keV (e.g. White & Holt 1982; Watanabe 
et al. 1994). The observed EFe spectrum in the high/soft 
state peaks at a few keV and is highly variable, whereas 
in the low/hard state the peak is at ~20 keV and there is 



less variability. At the other end of the spectrum, Cyg X- 
3 is the strongest radio source among the X-ray binaries 
and shows huge radio outbursts associated with relativis- 
tic jets. The radio activity and the production of the jets 
are closely linked with the X-ray emission and the differ- 
ent X-ray states (McCoUough et al. 1999; Choudhury et 
al. 2002). 

It is believed that Cyg X-3 shows so many unique prop- 
erties because it is in a very short-lived transitional phase 
of binary evolution. It is therefore possible that the sys- 
tem represents a stage of evolution that will be the fate 
of many of the massive X-ray binaries observed today. If 
this is the case, an understanding of this peculiar source 
could prove to have broad implications for high-mass X- 
ray binaries in general. 



2. OBSERVATIONS AND DATA ANALYSIS 



Cyg X-3 has been observed on several occasions with all 
three X/y-ray instruments aboard INTEGRAL. The most 
extensive observations were carried out during the INTE- 
GRAL performance verification phase observations of the 
Cygnus region in December 2002. (The data for Rev. 23 
are described in Vilhu et al. 2003, hereafter V03.) Cyg X- 
3 has also been in the FOV in several of the Galactic Plane 
Scans (GPS) carried out during 2003. Due to the large 
variability of the source (see Fig. 1), a simple co-adding 
of all available observations will not represent a physical 
spectrum of Cyg X-3. It is true, however, that the spectral 
state of this source can be roughly determined by the flux 
level in 1.5-12 keV X-rays as measured by the All Sky 
Monitor (ASM) on board RXTE. Here, we have therefore 
grouped the observations based on the RXTE/ ASM count 
rate and combined the selected science windows to build 
two co-added spectra of Cyg X-3, representing different 
spectral states. The dates for the observations are listed 
in Table 1 and also indicated in Fig. 1 together with the 
corresponding flux levels as seen by RXTE/ ASM. 
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Figure 1. Left panel: The RXTE/ASM lightcurve since launch ofRXTE in January 1996 until March 2004. The launch of 
INTEGRAL is shown by a vertical line. Right panel: An enlarged portion of the same lightcurve representing the lifetime 
of INTEGRAL. The times of the INTEGRAL observations are shown as arrows. Observations indicated with open arrows 
were used to construct the 'high ' state spectrum and observations indicated by filled arrows to construct the 'low' state ' 
spectrum. It is clear from a comparison with the long-term RXTE/ASM lightcurve that the source has not been in a very 
high state since the launch of INTEGRAL. 



Table 1. Data from the listed revolutions were used to 
create the 'high' state spectrum. The dates represent the 
start date of each revolution. The same dates are indi- 
cated by open arrows in Fig. 1. Effective exposure times 
in ksec are given for each instrument per revolution and 
in total. 



Table 2. Data from the listed revolutions were used to 
create the 'low' state spectrum. The dates represent the 
start date of each revolution. The same dates are indi- 
cated by filled arrows in Fig. 1. Effective exposure times 
in ksec are given for each instrument per revolution and 
in total. 



Rev. 


start date 


JEM-X 


ISGRI 


SPI 


19 


2002 - 12 - 09 


36 


36 




20 


2002 - 12 - 12 


17 


17 


133 


21 


2002 - 12 - 15 






180 


22 


2002- 12- 18 






161 


23 


2002 - 12 - 21 


82 


82 


169 


54 


2003 - 03 - 24 


2 


2 


9 


Total 


exposure time 


137 


137 


652 




2.1. 


JEM-X 







Rev. 


start date 


JEM-X 


ISGRI 


SPI 


59 


2003 - 04 - 08 


4 


2 




62 


2003 - 04 - 17 






11 


67 


2003 - 05 - 01 






12 


70 


2003 - 05 - 10 


2 


2 


11 


74 


2003 - 05 - 22 


4 


4 


11 


79 


2003 - 06 - 06 


4 


4 


12 


82 


2003 - 06 - 15 


2 


2 


11 


87 


2003 - 06 - 30 


2 


2 


14 


142 


2003 - 12 - 12 


2 


2 


11 


145 


2003 - 12 - 21 


2 


2 


11 


Total 


exposure time : 


22 


20 


104 



The observations were performed with the second (JEM- 
X2) of the two identical X-ray monitors. The JEM-X and 
ISGRI data were analysed using ISDC's (Courvoisier et 
al. 2003) Offline Standard Analysis (OSA) software ver- 
sion 3.0. The data were taken from 54 pointings from 
the high state and 1 1 pointings in the low state with 
the source always in the partially coded field-of-view, 
within 5° from the FOV center. Both Cyg X-1 and SAX 
J2 103. 5+4545 (an X-ray pulsar, discovered by Hulleman, 
in 't Zand & Heise 1998) are >10° from the FOV cen- 
ter and hence the data are not contaminated by those 
sources. Source spectra were extracted individually per 
pointing and then added together into an average 'high' 
and 'low' spectrum respectively. The spectra were im- 
plicitly background-subtracted by a deconvolution algo- 



rithm assuming a spatially flat background. We used the 
energy range of 3.0-20 keV for spectral fitting. System- 
atic errors were added as 20% for channels 1-58, 10% 
for channels 59-96 and 2% for channels 96-255. 



2.2. IBIS/ISGRI 

For ISGRI, we have selected 54 INTEGRAL pointings 
for the high state and 10 for the low state with the source 
always in the totally coded field-of-view. This way we 
have simultaneous ISGRI and JEM-X data and we do 
not use the ISGRI partially coded field-of-view for which 
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the ISGRI response matrix is known to be less accurate. 
The two different spectral-state set of pointings have been 
analysed separately. For each set we have produced a 
mosaicked image (weighted combination of single point- 
ing images) in 10 energy bands. We built the high and 
soft spectrum extracting the counts/sec from the bright- 
est pixel in each energy band. This method of extracting 
the spectra from the mosaics instead of adding individual 
spectra significantly reduces the level of noise and gives 
high-quahty spectra even for low count rates. For the fit- 
ting of the two broad-band spectra, we used the energy 
range 20-150 keV. A 5% systematic error was added to 
all channels. 



2.3. SPI 

SPI observations within 15° from the FOV center, from 
times with similar flux levels at 1.5-12 keV as for the 
JEM-X and ISGRI data, were selected. The SPI data 
were analysed using OSA version 3.0, except for some 
updates due to known problems in the binning software. 
For image reconstruction and spectral extraction we ap- 
plied version 6.0 of the SPI Iterative Removal Of Sources 
program (SPIROS; Skinner & ConneU 2003). This al- 
lowed us to apply a new background model based on the 
mean count modulation of the detector array. In order 
to get precise fluxes, the source positions of the known 
sources in the field of view have been fixed to their cat- 
alogue values. In addition, we allowed SPIROS to apply 
time dependent normalisation to the source fluxes of the 
three brightest and variable sources (Cygnus X- 1 , Cygnus 
X-3, and EXO 2030H-375). Thus the SPI spectrum of 
Cygnus X-3 presented here is an average over 10 ksec 
time bins. For spectral extraction, 25 logarithmic bins in 
the 20-100 keV energy range, and 5 logarithmic bins in 
100-2000 keV have been used. The energy range 20- 
200 keV was used for the spectral fitting due to strong 
background features at higher energies. The instrumental 
response function has been derived from on-the-ground 
calibration (Sturner et al. 2003) and then corrected based 
on the Crab calibration observation. A 5% systematic er- 
ror was added to all channels. 

The resulting two combined JEM-X, ISGRI and SPI 

broad-band spectra were fitted using the XSPEC package 
(Arnaud 1996). The ISGRI and SPI spectra have been 
renormalized to the JEM-X data. The ISGRI and JEM-X 
levels showed good agreement while the SPI normahza- 
tion is a factor of ~2 higher. The obtained 'high' spectra 
correspond to the 1.5-12 keV source flux >200 mCrab, 
and the 'low' spectra, to <130 mCrab. Note that neither 
of these states represent the extreme high/soft or low/hard 
states displayed by this source on a longer timescale (see 
§4 below). 



3 . BROAD-BAND SPECTRAL MODELLING 



In spite of a large number of X-ray observations, the spec- 
tra of Cyg X-3 have mostly been interpreted in terms of 
phenomenological models (e.g. Nakamura et al. 1993), 



such as power law and bremsstrahlung. The first physi- 
cal interpretation of the intrinsic broad-band spectrum of 
Cyg X-3, based on simultaneous INTEGRAL (Rev. 23) 
and /^ZrS/PCA-HEXTE observations in 2002 Dec. 22- 
23, was given by V03. Their model assumes a hot Comp- 
tonizing plasma in the vicinity of an accretion disc, which 
both provides soft seed photons for the Comptonization 
and reflects the hard photons emitted by the hot plasma. 
The specific model, which we also use in this work, is 
eqpair (Coppi 1992, 1999). 

In general, the electron distribution in the hot plasma can 
be purely thermal or hybrid, i.e., Maxwellian at low ener- 
gies and non-thermal at high energies, if an acceleration 
process is present. This distribution, including the elec- 
tron temperature, T, is calculated self-consistently from 
the assumed form of the acceleration (if present) and from 
the luminosities corresponding to the plasma heating rate, 
Lh, and to the seed photons irradiating the cloud, Lg. The 
plasma optical depth, t, includes a contribution from e* 
pairs. The importance of pairs depends on the ratio of 
the luminosity to the characteristic size, r, which is usu- 
ally expressed in dimensionless form as the compactness 
parameter, ( = LcTTlirm^c'), where crx is the Thomson 
cross section and me is the electron mass. 

In the fitting process, we have found that our both broad- 
band INTEGRAL spectra are compatible with the hot 
plasma being completely thermal, and kT « 511 keV. 
Then, the e* pair production is negligible, and the ab- 
solute value of the compactness is only weakly impor- 
tant. Accordingly, we assume a constant 4 = 10 (as 
in V03). We stress, however, that our data, limited to 
energies <200-300 keV, rather weakly constrain possi- 
ble non-thermal component of the electron distribution 
and the presence of pairs. Non-thermal processes cer- 
tainly exist in at least some states of Cyg X-3 (Szostek & 
Zdziarski 2004; see §4 below). 

Our data, covering only photon energies >3 keV, rather 
poorly constrain the form of the seed photons. We there- 
fore use here the same simple assumption as in V03, i.e, 
emission of an accretion disc with the maximum temper- 
ature fixed in the fit. The data also do not allow to con- 
strain the structure of the X-ray absorber and the likely 
presence of additional spectral components in soft X-rays 
(e.g. Molnar & Mauche 1986; Nakamura et al. 1993). 
This absorption and emission is due to the dense wind 
from the companion Wolf-Rayet star. These processes 
(which we do not treat in detail here) may strongly affect 
the spectra up to as much as ~20 keV. 

For the purpose of modelling the broad-band spectra, we 
use the same phenomenological model of the absorber 
as V03. It assumes that one absorbing medium with the 
column density, A'h,o> fully covers the source, and an- 
other medium with the column, A^h,i, covers a fraction, 
/i, of the source. We also include Compton reflection 
(Magdziarz & Zdziarski 1995), parametrized by an effec- 
tive solid angle subtended by the reflector as seen from 
the hot plasma, Q, and assuming an inclination of 60°. 
We also include an Fe Ka fluorescent line, which we 
model as a Gaussian (with the width of cr-Ka, the pho- 
ton flux of Fy^a, and the peak energy at E^„). We al- 
low the reflecting medium (at the assumed temperature 
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Table 3. Model parameters'^ for the two broad-band INTEGRAL spectra of Cyg X-3 in the 'high' and 'low' state. 



Data 


10^2 cm-- 


10^2 cm-2 


/l 


4/fs 


T 


keV 


n/2;r 


erg cms"' 


keV 


-'' -1 
cm - s 


ergcm"2 §-1 


x-ly 


High 
Low 


20.7 ±1.4 
8.6±7.3 


570±50 
32±21 


0.14±0.11 
0.75±0.20 


0.21 ±0.03 
0.22±0.06 


0.21 ±0.04 
0.15±0.13 


73 
102 


0.6±0.4 
2.0±1.2 


7000+3000 
'"""-7000 

3000 ±1000 


6.61±0.12 
6.87±0.07 


0.0085 
0.0146 


5.2x10-' 
3.8x10-' 


129/172 
129/147 



" The uncertainties are the linear uncertainties as given by XSPEC; the blackbody temperature was fixed at 0.38 keV. 
Calculated from the energy balance, i.e., not a free fit parameter. 
Assumed < iC in the fits. 

The line photon flux corrected for the Galactic N\{ only. 

The bolometric flux of the absorbed model spectrum and normalized to the JEM-X/ISGRI data. 



T i 1 — I — I — r-r-| 1 1 1 1 — I — I — r-r-| 1 r 




E [keV] 

Figure 2. Deconvoled spectra of Cyg X-3. The green, red and black data points correspond to the JEM-X, ISGRI and SPI 
spectra, respectively, in the high state. The heavy blue, cyan and magenta data points are for the JEM-X, ISGRI and SPI, 
respectively, in the low state. The model spectra are shown in the green and blue curves. The ISGRI and SPI spectra have 
been renormalized to the JEM-X data. The proportions of this figure correspond to equal length per decade on each axis. 
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of 10^ K) to be ionized, using the ionization calculations 
of Done et al. (1992). The ionizing parameter is defined 
as ^ s AnFionIn, where Fion is the ionizing flux and n 
is the reflector density. Given the simplified treatment of 
the ionized reflection of Done et al. (1992), we assume 
^ < lOl 

The obtained model parameters are given in Table 3, and 
the deconvolved data and the model spectra are shown in 
Fig. 2. We see that the two spectra differ strongly at low 
energies, but their high-energy tails, at energies >20 keV, 
are almost identical, with only a slight difference in the 
normalization. 



4. DISCUSSION 



Somewhat surprisingly, the two average spectra obtained 
by us are very similar at energies above 20 keV, even if 
there is a clear difference between the shape below 20 
keV. We can compare them to the range of the spectra ob- 
served so far by the pointed instruments on board RXTE. 
A study of the different spectral states observed by the 
PCA and HEXTE was performed by Szostek & Zdziarski 
(2004), who used the same model as that used by us and 
by V03. Fig. 3 shows four out of five average spectra 
of Szostek & Zdziarski (2004), with the second hardest 
spectrum (group #2) omitted for clarity. Those spectra 
have also been shown by Zdziarski & Gierlihski (2004) 
to be similar to those of the canonical states of black-hole 
binaries (notwithstanding the peculiar form of absorption 
in CygX-3). 

In Fig. 3, we also show the observed 'low' and 'high' 
INTEGRAL spectra. We find them to be relatively similar 
in shape to the hard (blue) and intermediate (magenta) 
average RXTE spectra. However, the normalization of 
the INTEGRAL spectra is lower by a factor of ~2. A part 
of this normalization difference is likely to be due to the 
different absolute calibration of the JEM-X and ISGRI 
vs. that of the PCA. 

From Fig. 3, we can clearly conclude that with INTE- 
GRAL, we have as yet sampled neither the hardest nor the 
softest spectra of Cyg X-3. Of particular interest here is 
the ultrasoft state (red spectrum in Fig. 3) with a distinct 
hard tail. The presence of such a high energy tail with- 
out a cutoff below ~200 keV requires a nonthermal elec- 
tron distribution and an acceleration mechanism operat- 
ing. Due to the fact that the HEXTE detector can mea- 
sure it only up to ~ 100 keV, the detailed form of this hard 
tail in the high/ultrasoft state has remained largely unex- 
plored, and it would be of great importance to observe 
it with INTEGRAL. The presence of that tail is clearly 
accompanied by the very highest level of the flux at 1 .5- 
12 keV, which is also monitored by the RXTE/ ASM. All 
ultrasoft spectra corresponding to the average shown in 
Fig. 3 were observed when the ASM flux level was >27 
counts/sec, which corresponds to >360 mCrab. In com- 
parison, see the right panel of Fig. 1, the highest flux 
levels observed during the lifetime of INTEGRAL have 
only occasionally reached above 20 counts/sec and IN- 
TEGRAL has thus yet not been able to observe a true 
high/ultrasoft state of Cyg X-3. 
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Figure 3. The range of spectra observed so far by 
RXTE/PCA/HEXTE compared to the present INTEGRAL 
spectra. The red, magenta, green and blue spectra 
represent four out of five groups obtained by Szostek 
& Zdziarski (2004) by averaging 42 RXTE observa- 
tions from 1996-2000 (normalized to the PCA). The low 
and high INTEGRAL spectral models are plotted in the 
dashed and dotted cyan curves, respectively. 



5. CONCLUSIONS 



We have combined multiple JEM-X, ISGRI and SPI 
spectra from available public and GPS observations to 
create two broad-band spectra of Cyg X-3, correspond- 
ing to high and low levels of the 1.5-12 keV flux. The 
resulting 'high' and 'low' spectra resemble what can be 
referred to as the intermediate and the low/hard state of 
this source, but are similar in shape above 20 keV. We 
have fitted the intrinsic spectra with thermal Comptoniza- 
tion models including Compton reflection, which are (see 
Zdziarski & Gierlinski 2004) similar to those found for 
black-hole binaries. We point out that INTEGRAL has 
yet to observe Cyg X-3 in the high/ultrasoft state, which 
would provide a very interesting opportunity to study the 
hard non-thermal tail associated with this state. 
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